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ABSTRACT 


The  estimated  weight  advantage  of  designing 
a  superconducting  homopolar  acyclic  machine  with 
a  saturated  magnetic  shield  is  presented,  and  the 
possible  power  or  flux  loss  due  to  the  saturation 
of  the  iron  is  reported.  A  method  of  estimating 
the  magnetic  field  at  15.24  cm  (6  in)  from  the 
saturated  magnetic  shield  is  shown  in  addition  to 
indicating  how  the  fields  at  large  distances  vary 
for  hexapole  and  quadrupole  designs. 
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INTRODUCTION 


Superconductive  acyclic,  or  homopolar,  machines  are  being  developed 
by  the  Navy  for  ship  propulsion  application.  The  machine  contains  three 
major  elements;  a  superconducting  coil  to  generate  mangetic  flux,  an 
armature  circuit  which  converts  electrical  power  to  mechanical  power  in 
the  presence  of  a  magnetic  field,  and  a  ferromagnetic  shield  to  shape 
and  contain  the  flux.  A  machine  may  include  a  single  solenoid  (dipole) 
magnet  winding,  as  illustrated  conceptually  in  Figure  1,  or  employ  two 
or  three  solenoids  (i.e.  quadrupole  or  hexapole  windings).  In  Figure  1, 
the  superconducting  magnet  within  the  helium  vessel,  or  dewar,  is  the 
innermost  machine  element.  The  intense  magnetic  flux  generated  in  the 
bore  of  the  solenoid  is  attracted  by  the  ferromagnetic  shield  forcing 
virtually  all  of  the  magnetic  flux  to  radially  transverse  the  rotor 
twice.  When  current  is  passed  through  brushes  and  axially  down  the 
copper  rotor  conductors,  the  resulting  Lorentz  interaction  provide  motor 
action. 

In  the  design  of  superconducting  homopolar  machines,  the  practice 
has  been  to  design  the  magnetic  shield  in  the  unsaturated  state,  15-20 
kilogauss  (KG).  That  is,  to  design  the  magnetic  shield  so  that  the  flux 
leakage,  close  to  the  machine  surface  (typically  15.24  cm  (6  in)  from 
the  shield),  doesn't  exceed  a  flux  density  of  100  gauss.  The  low  flux 
leakage  was  established  so  that  the  magnetic  field  generated  by  the 
superconducting  magnet  would  not  effect  the  operation  of  near-by¬ 
equipment. 

However,  it  was  found  that  relaxation  of  the  flux  leakage 
requirement  could  result  in  a  very  large  machine  weight  reduction  due  to 
the  mass  density  and  volume  of  the  ferromagnetic  shield.  This  study  was 
undertaken  to  quantify  this  potential  machine  weight  reduction. 

In  order  to  establish  the  weight  benefit  of  relaxing  the  flux 
leakage  requirements,  an  analysis  was  conducted  and  the  machine  weight 
advantage  evaluated  for  a  specific  machine  design  and  power  rating  (i.e. 
a  1.27  weber  hexapole  design).  In  addition,  a  scaling  method  was 
established  to  extend  the  weight  reduction  estimates  for  a  specific  flux 
leakage  to  machines  of  different  power  ratings.  The  field  values  for 
the  hexapole  and  quadrupole  machine  at  large  distances  are  shown  to  vary 
as  1/R3  and  1/R  ,  respectively. 
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Shaped  Field  Machine  Concept 


The  present  analysis  will  cover  the  following  areas: 

o  The  weight  advantage  available  due  to  the  allowable  flux 
leakage. 

o  An  evaluation  of  the  power  loss  due  to  the  flux  leakage. 

o  An  evaluation  of  a  linear  scaling  equation  for  the  flux-leakage 
near  the  machine  (i.e.  15.24  cm  (6  in)  from  the  outside 
diameter,  OP). 

o  The  scaling  of  fields  at  large  distances  from  different  machine 
designs. 

o  The  comparison  of  the  flux  leakage  for  a  shielded  and  an 
unshielded  (self-shielded)  machine. 


ANALYSIS 

In  this  study,  two  computer  programs  were  used  for  calculating  the 
magnetic  fields  of  these  machines.  The  first  is  a  2-dimensional 
triangular  irregular  mesh  (TRIM)  computer  program  with  the  capabilities 
of  calculating  magnetic  fields  in  the  presence  of  iron.  TRIM  was 
developed  at  Arqonne  National  Laboratory.  The  second  computer  program 
(IMAGF)  uses  the  Biot-Savart  Law  to  calculate  the  magnetic  field  at 
specified  "mesh  points"  due  to  a  number  of  "coil  regions",  and  assumes 
no  iron  regions. 

The  configuration  of  the  hexapole  machine  design  used,  and  the 
region  of  interest  for  the  flux  leakage,  is  shown  in  Figure  2.  An 
average  flux  density  of  15  KG  in  the  shield  of  a  hexapole  design  was 
used  as  the  base  case  (TRIM  run  SMR  97B).  Using  TRIM,  a  percentage  of 
the  magnetic  shield  was  removed  and  the  flux  leakage  was  noted  at 
several  locations.  This  information  is  presented  in  Table  1  and  the 
flux  plots  for  several  of  these  cases  are  shown  in  Figures  3.  The 
graphical  representation  of  this  data  is  shown  on  Figures  4  and  5, 
showing  the  potential  weight  reduction  due  to  the  increase  in  flux 
leakage.  In  addition.  Figures  4  and  5  indicate  that  designing  a  shield 
at  saturation  (20  KG)  provides  a  20 %  decrease  in  the  shield  weight  with 
very  little  flux  leakage.  This  becomes  significant,  when  it  is  realized 
that  the  weight  of  the  shield  in  the  larger  machines,  40,000  horsepower 
(HP),  contribute  about  50-80%  of  the  total  weight  of  the  machines.  In 
Figure  4,  field  values  at  the  larger  radii  were  calculated  based  on  the 
assumption  that  they  varied  as  1/R  ,  which  will  be  shown  later  to  be 
correct  for  a  hexapole  design. 
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Weight  Reduction  of  the  Magnetic  Shield  Vs.  Flux  Leakage 
at  Various  Locations  for  1.27  Weber  Hexapole  Design 


The  removal  of  a  portion  of  the  magnetic  shield  has  been  shown  to 
have  quite  a  weight  benefit.  However,  it  also  has  a  disadvantage  of 
increasing  the  reluctance  of  the  magnetic  circuit  with  the  consequent 
decrease  flux  and  machine  power,  as  indicated  in  Figure  5  and  Figure  6. 
In  fact,  the  removal  of  the  magnetic  shield  reduces  the  effective  flux 
by  approximately  25%  at  constant  field  excitation.  This  decrease  in 
power  or  flux  can  be  easily  compensated  by  increasing  the  magnet 
excitation  (i.e.  amp-turns)  slightly,  resulting  in  less  than  a  2% 
increase  in  the  machine  diameter. 

In  order  to  evaluate  the  approximated  flux  leakage  density  at  a 
specific  radial  distance  (15.24  cm  (6  in)  from  outside  diameter,  OD,  of 
the  machine),  a  linear  scaling  technique  was  structured.  The  technique 
requires  knowledge  of  the  field,  flux  and  geometry  of  a  reference  coil 
(coil  1)  to  estimate  the  field  characteristic  of  a  second  coil  (coil 
2).  This  scaling  equation  for  a  hexapole  or  quadupole  machine  is 
defined  as: 


b2  =  B1  ^e2/0el)  {trl2  "  roi2Mr22  -  r02])  (1) 

The  scaling  equation  used  in  this  study  was  derived  based  on  a  1.27 
weber  hexapole  design  (SMR  97H)  and  then  used  to  predict  the  flux 
leakage  at  higher  flux  level  hexapole  machines  (i.e.  up  to  5.0  webers) 
and  the  predictions  of  hexapole  machines  with  varying  amount  of  magnetic 
shielding.  Since  the  scaling  equation  used  in  this  study  was  linear, 
the  non-linear  effects  become  apparent  when  the  results  were  compared  to 
the  TRIM  output.  This  non-linear  effect  is  illustrated  in  Figure  7. 

The  linear  scaling  technique  developed,  in  equation  1,  appear  to  over 
predict  at  higher  flux  levels  and  under  prediction  for  reductions  in  the 
shielding.  A  more  accurate  prediction  could  be  made  by  adding  to  the 
scaling  equation,  a  function  that  would  empirically  account  for  some  of 
these  non-linear  effects.  The  linear  correction  factor  (Fc2)  for 
scaling  up  (i.e.  coil  2  is  larger  than  coil  1)  in  flux  from  the  hexapole 
design  SMR  97H  is  defined  as: 

Fcl  =  -.05404  0e  +  1.0669 

The  correction  factor  (F_2)  for  scaling  down  in  shielding  from  the 
hexapole  design  SMR  97H  Vs  defined  as: 

Fc2  =  7.5198  x  10"3  E  +  .70498 

These  linear  correction  factors  were  than  factored  into  the  scaling 
equation  and  using  the  5.0  weber  hexapole  design,  the  flux  density 
leakage  at  15.24  cm  (6  in)  was  estimated.  Results  showed  a  maximum 
error  of  less  than  9%. 
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5.0  WEBER  HEXAPOLE  DESIGN 
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%  WEIGHT  REDUCTION  OF  SHIELD 
Figure  6 


WEIGHT  REDUCTION  OF  SHIELD 

Figure  7 


The  computer  solution  for  the  hexapole  magnet  system  was  the 
compared  to  the  functions  1/R  ,  1/R  ,  and  1/R  ,  and  the  results  are 
shown  on  Figure  8.  It  was  found  that  no  good  comparison  resulted  from 
the  curves  of  the  magnetic  field  vs.  radial  distance.  However,  it  was 
found  that  at  these  small  distances  from  the  machine,  the  field  vaules 
seemed  to  all  fall-off  faster  than  1/R3.  As  a  result,  the  assumption 
that  the  leakage  flux  for  a  hexapole  magnet  deisgn  will  vary,  at  the 
very  worst,  as  1/R3.  Using  the  computer  program  IMAGE  for  the  case  of 
no  magnetic  shield  it  was  found  that  at  large  radii  the  field  varies 
very  closely  as  1/R3  for  a  hexapole  design  (see  Figure  9).  It  was  also 
noted  that  at  large  radii,  the  fields  of  a  solenoid  varies  as  1/R3  (see 
appendix),  and  the  fields  of  a  quadrupole  design  varies  as  1/R**  (see 
Figure  10).  Thus,  showing  that  a  quadropole  machine  design  has  an 
advantage  over  a  solenoid  or  a  hexapole  design  in  decreasing  its  flux 
leakage  at  large  radii. 

A  comparison  was  made  between  a  shielded  and  an  unshielded  hexapole 
and  quadrupole  machines  in  order  to  determine  how  the  flux  leakage 
varied  with  machine  0.0.  The  results  for  a  hexapole  deisgn  is  indicated 
on  Figure  11  and  the  quadrupole  on  Figure  12.  Figure  11  shows  a  dip  in 
the  field  values  for  a  hexapole  design,  this  is  believed  to  be  caused  by 
the  interaction  of  the  end  coil  field  with  the  center  coil  field.  In 
addition,  it  may  be  noted  that  a  much  greater  decrease  in  the  field 
value  occurs  for  the  shielded  hexapole  machine  than  the  unshielded 
case.  This  large  decrease  in  field  value  is  also  noted  in  Figure  12  and 
represents  the  location  of  the  outer  surface  of  the  magnetic  shield. 

CONCLUSIONS 

The  following  conclusions  were  drawn  from  this  study: 

o  A  very  large  weight  advantage  can  be  realised  if  the  restriction 
on  the  flux  leakage  (100  gauss  at  15.24  cm  (6  in)  from  the 
machine)  is  reduced. 

o  The  power  loss  due  to  the  removal  of  the  shield  is  very  small 
and  can  be  made  up  by  increasing  the  diameter  of  the  machine  by 
less  than  27,. 

o  At  large  radii  the  field  values  varies  as: 

1/R3  for  a  hexapole  design 
1/R4  for  a  quadrupole  design 
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FIELD  VALUE  (GAUSS) 


?  RADIAL  DISTANCE  (CM) 
Figure  8 


Comparison  of  the  Computer  Program  TRIM  and  the  Scaling  Technique 
for  a  Hexapole  Magnet  System 
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FiELD  VALUE  (GAUSS) 


RADIAL  DISTANCE  (cm) 

Figure  9 

Comparison  of  the  Computer  Program  IMAGE  and  the  Scaling  Technique  for  a  Hexapole 

Magnet  System  (No  Iron  Shielding) 
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RADIAL  DISTANCE  (NO.  OF  OD's) 
Figure  11 


Comparison  of  the  Magnetic  Fields  of  a 
Shielded  and  Unshielded  Hexapole  Machine 
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ABSOLUTE  MAGNETIC  FIELD  (GAUSS) 


Comparison  of  the  Magnetic  Fields  of  a  Shielded  and  Unshielded 

Quadrapole  Machine 
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APPENDIX  A 


MAGNETIC  FIELDS  FOR  A  SOLENOID,  AT 
LARGE  DISTANCES 
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APPENDIX  A 

MAGNETEIC  FIELDS  FOR  A  SOLENID.  AT 

- LARGE  DISTANCE'S- - - 

This  derivation  indicates  that  fields  from  a  solenoid  having  no 
iron  varies  as  1/R  .  The  basic  deviation  was  obtained  from  The 
Introduction  to  Electromagnetic  Fields  and  Waves,  by  Charles  H.  Holt, 
pages  366-358. 

The  solenoid  and  the  coordinate  systems  are  shown  in  Figure  A-l: 


z 


where  vector  potential 

J  *  current  density  vector 
u*  permeability  of  the  media 
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A-2 


In  the  RECTANGULAR  COORDINATE  SYSTEM,  this  equation  has  the  same  form  as 
Poisson's  equation.  The  integral  solution  of  Equation  A-l  becomes: 


A  =  dv 

4tt  v  r 


(A— 2 ) 


In  terns  of  the  steady  current  I  of  a  filamentary  conductor,  the  vector 
potential  becomes: 


A  = 


57 


Id* 


c  r 


(A— 3 ) 


It  is  noted  from  symmetry  considerations  that  the  magnitude  of  the 
vector  potential  (^)  is  independent  of  0.  Therefore,  0  can  be  set 
at  tt/2  resulting  in  only  A  X-directional  component  of  the  vector 
potential  =  -A„).  In  addition,  because  of  the  orientation  of  the 
coil  and  the  coordinate  system  (spherical)  used,  the  vector  potential 
has  only  the  single  spherical  component  A0. 


Therefore,  at  the  Point  P, 


-A 


x 


Thus,  Equation  (A-3)  can  be  written  in  the  rectangular  coordinate  system 
in  the  following  manner: 


Ax 


JL  f  I 

4n  c  r 


dx 


(A— 4 ) 


Z 
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Equation  (A-4)  can  then  be  put  in  a  sperical  coordinate  system  in  the 
following  manner:  (See  Figure  A-2  above) 


dS.  ^ado 

. _ .  -dx 

cos  (90-0)  *  diT 

di,  =  -dx/cos  (90-0) 

Using  the  identity  sin  0  =  cos  (90-0)  results  in  the  following: 
d£  =  -dx/sin  0 
ad0  =  dx/sin  0 


dx  =  -a  sin  0  d  0  (A-5) 

Substituting  equation  (A-5)  into  equaton  (A-4)  to  give  the  vector 
potential  gives: 


A0  =  Ax 

A0  =  -  JL  /  I  (-a  sin  0)  d  0  (A-6) 

4ir  r 

A0  =  yla  /2tt  sin  0  d  0 
4tt  o  r 

Equation  (A-6)  has  the  variable  (r)  which  is  a  function  of  0.  Thus,  r 
can  be  rewritten  in  the  following  manner  (see  Figure  A-l): 

r2  =  (X!  -  x)2  +  (y  -  yj)2  +  Z2 

r2  =  Xj2  _  2xjx  +  x2  +  y2  -  2yyl  +  y  2  +  z2 

It  is  noted  from  Figure  A-l  that  a  and  rQ  can  be  written  as: 

a2  =  xj2  +  yi2 

r02  3  y2  +  Z2 

Substituting  the  relations  for  a2  and  r2  in  the  equation  for  r2  gives: 

r2  =  a2  +  r02  -  2rQ  a  sin0  sin0 
and 

r  »/a2  +  r02  -  2r0  a$1n  6  sin  0 

substituting  the  above  relation  for  r  Into  equation  (A-6)  results  in 
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sin  0  d  0 


2n 


A0 


yla  j-  - 

4lr  0  /  a2  +  r  2  -  2rn  a  sin  e  sin  0 


(A-7) 


This  equation  is  very  hard  to  solve  because  of  the  use  of  the  elliptical 
integrals.  However,  if  the  following  is  assumed 


a  «  rr 


the  solution  becomes  greatly  simplified.  Then  using  this  assumption  in 
equation  (A-7)  gives  the  following 


A0 


yla  . 

4tt 


2  u 


sin  0  d  0 


0  /  1  -  (23/r0 )  sin  e  sin  0 


(A-8) 


Then  noting  that  the  second  term  in  the  radical  is  very  small  compared 
with  unity,  the  denominator  can  be  approximated  by  using  a  binominal 
expansion  and  neglecting  all  except  the  first  two  terms.  Thus, 

(1  -  (2a/r0)-*^  sin  e  sin  0x^1+-^-  sin  e  sin  0 

When  substituted  into  equation  (A-8) 

A0  =  ■  f  (sin  0)  (1  +  ~  sin  e  sin  0)  d  0 


4?r  rQ  o 

Integrating  equation  (A-9)  gives 


(A-9) 


A  0 


ylS  sin  e 
4tt  r^ 


where  S  =  ira* 


The  magnetic  field  (B)  can  be  evaluated  by  taking  the  curl  of  A0 
B  =  V  X  A0 


n _ X—  r_l  (A.  sin  e)  -  i 

°  "  r  sin  0  ^30  to  30  J 


+  -°  [  -J— 

r  L  sin  e 


3A 
_ r 

30 


3r 

o 


( r A0 ) 3 


(A-10) 
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